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Abstract 

We study the phases of strongly correlated electron systems in two spatial dimensions in the 
framework of AdS4/CFT3 correspondence. The AdS (gravity) model consists of a Dirac fermion 
coupled to electromagnetic field and gravity. To classify the ground states of strongly correlated 
electrons on the CFT side and to construct the full phase diagram of the system, we construct a 
quantum many-body model of bulk fermion dynamics, based on the WKB approximation to the 
. Dirac equation. At low temperatures, we find a quantum corrected approximation to the electron 

I star where the edge is resolved in terms of wavefunctions extended fully through AdS. At high 

temperatures, the system exhibits a first order thermal phase transition to a charged AdS-RN 
black hole in the bulk and the emergence of local quantum criticality on the CFT side. This 
change from the third order transition experienced by the semi-classical electron star restores the 
intuition that the transition between the critical AdS-RN liquid and the finite density Fermi system 
is of van der Waals liquid-gas type. 
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I. INTRODUCTION 



The problem of fermionic quantum criticality has proven hard enough for condensed 
matter physics to keep seeking new angles of attack. The main problem we face is that 
the energy scales vary by orders of magnitude between different phases. The macroscopic, 
measurable quantities emerge as a result of complex collective phenomena and are difficult 
to relate to the microscopic parameters of the system. An illustrative example present the 
heavy fermion materials [sl which still behave as Fermi liquids but with vastly (sometimes 
hundredfold) renormalized effective masses. On the other hand, the strange metal phase 
of cuprate-based superconducting materials j4|, while remarkably stable over a range of 
doping concentrations, shows distinctly non-Fermi liquid behavior. The condensed matter 
problems listed all converge toward a single main question in field-theoretical language. It 
is the classification of ground states of interacting fermions at finite density. 

In this paper we attempt to understand these ground states in the framework of 
AdS/CFT, the duality between the strongly coupled field theories in d dimensions and 
a string configuration in d + 1 dimension. Holography (AdS / CFT correspondence) 
has become a well-established treatment of strongly correlated electrons by now, but it still 
has its perplexities and shortcomings. Since the existence of holographic duals to Fermi 
surfaces has been shown in , the next logical step is to achieve the understanding of 
the phase diagram: what are the stable phases of matter as predicted by holography, how 
do they transform into each other and, ultimately, can we make predictions on quantum 
critical behavior of real- world materials based on AdS/CFT. 

The classification of ground states now translates into the following question: classify 
the stable asymptotically AdS geometries with charged fermionic matter in a black hole 
background. Most of the work done so far on AdS/CFT for strongly interacting fermions 
relies on bottom-up toy gravity models and does not employ a top-down string action. We 
stay with the same reasoning and so will work with Einstein gravity in 3 -|- 1 dimensions. 
We note, however, that top-down constructions of holographic fermions exist gI, 

So far three distinct models aiming at capturing the stable phases of holographic fermionic 
matter have appeared: the electron star 10|, Dirac hair 13| and a confined Fermi liquid 
model [l^. The electron star is essentially a charged fermion rewriting of the well-known 
Oppenheimer-Volkov equations for a neutron star in AdS background. The bulk is thus 



2 



modeled as a semiclassical fluid. The mystery is its fleld theory duah it is a hierarchically 



ordered (inflnite) multiplet of fermionic liquids with stable quasiparticles [l8j. On the other 
end of the spectrum is Dirac hair, which reduces the bulk fermion matter to a single radial 
harmonic. As a consequence the fleld theory dual is a single Fermi liquid, however its 
gravitational consistency properties are not yet fully understood. In 19| we have shown 
that Dirac hair and electron star can be regarded as the extreme points of a continuum 
of models, dialing from deep quantum - a single radial mode — to a classical regime — a 
very large occupation number — in the bulk. They correspond to two extreme "phases" in 
the fleld theory phase diagram: a multiplet of a very large number of Fermi liquids and a 
single Fermi liquid. The third model, the conflned Fermi liquid [l^ introduces conflnement 
through modifying the bulk geometry and solves for quantum-mechanical wave functions 
adding them up to compute the full bulk density. This latter step is more general then the 
single-particle approach of jisl and it naturally extends the single harmonic Dirac hair state 
with a single Fermi surface to a state with multiple Fermi surfaces. Our main motivation is to 
construct a complementary model that extends from the other end — the semi-classical fluid 
— down to a state with a countable number (a few) Fermi surfaces. We aim for a system 
which is general enough to encompass the middle ground between extreme quantum and 
extreme classical regimes in the original deconflned setup. In the recent model of "quantum 
electron star" [l^] the same goal is set but the method used is different and is based on the 
deconflned limit of 1^ . 

In addition to simply improving the mathematical treatment of the bulk many-body 
fermion system, the guiding principle in our analysis will be to rest on the advantages and 
disadvantages of the current models. On the one hand, the Dirac hair is a fully quantum- 
mechanical model which shows its strength in particular near the boundary (the ultraviolet 
of the fleld theory) but becomes worse in the interior, i.e. close to the horizon (the infrared 
of the fleld theory) where density is high and the resulting state of matter cannot be well 
described by a single wavefunction. On the other hand, the electron star yields a very robust 
description of high-density matter in the interior but its sharp boundary at some radius Vc is 
clearly incompatible with a fully quantum description. It is thus obvious that the physically 
interesting model lies somewhere in-between the two approaches. 

We will use a simple WKB formalism to approximate the many-body Fermi system in 
the AdS bulk. This adds quantum corrections to the Thomas-Fermi (fluid) approximation 
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by taking into account finite level spacing. The occupation number itself acts as the control 
parameter of our approximation. The most notable feature, however, occurs in the transi- 
tion from the semiclassical approximation at infinite occupation number to finite occupation 
number. We find that the finite density quantum many body phases with fermionic quasi- 
particles at high enough temperatures always exhibit a first order transition into the zero 
density AdS-RN phase. Intuitively, this can be interpreted as a universal van der Waals 
liquid-gas transition. On the other hand in the semiclassical fluid limit underlying the elec- 
tron star, the transition was found to be continuous 
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34| . With this re-emergence of the 



first order nature of the thermal phase transition at the quantum level our results confirm the 
intuition that a density driven phase transition is always first order as also indicated by the 



Dirac hair approximation 



13| . One should take home that in the context of fermionic ques- 



tions in AdS/CFT quantum "1/N" corrections can be important and that the semiclassical 
fluid limit can be unreliable. 

The outline of the paper is as follows. In the Section II we describe the field content 
and geometry of our gravity setup, an Einstein-Maxwell-Dirac system in 3 + 1 dimension, 
and review the single-particle solution to the bulk Dirac equation. In Section III we start 
from that solution and apply the WKB approximation to derive the Dirac wave function 
of a many-particle state in the bulk. Afterwards we calculate density and pressure of the 
bulk fermions - the semiclassical estimate and the quantum corrections, thus arriving at the 
equation of state. Section IV contains the solution of the self-consistent set of equations 
for fermions, gauge field and the metric. There we also describe our numerical procedure. 
Section V is the core, where we analyze thermodynamics and spectra of the field theory side 
and identify different phases as a function of the three parameters of the system: chemical 
potential /i, fermion charge e and conformal dimension A. Section VI sums up the conclu- 
sions and offers some insight into possible broader consequences of our work and into future 
steps. 



II. HOLOGRAPHIC FERMIONS IN CHARGED BACKGROUND 

We wish to construct the gravity dual to a field theory at finite fermion density. We will 
specialize to 2 + 1-dimensional conformal systems of electron matter, dual to AdS4 gravities. 
We consider a Dirac fermion of charge e and mass m in an electrically charged gravitational 
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background with asymptotic AdS geometry. Adopting the AdS radius as the unit length, 
we can rescale the metric g^^ and the gauge field A^: 



(1) 



In these units, the action of the system is: 

1 



S 



—L'{R + 6) + -F' + L'C, 



T 



(2) 



where k is the gravitational coupling and F^i, = d^Ay — d^A^ is the field strength tensor. 
The fermionic Lagrangian is: 

1 



/ 



e^r^ ( + ^ouI^^Tbc - leLA^ ) - ml 



(3) 



where = i'^^T^, is the vierbein and is the spin connection. 

We shall be interested in asymptotically AdS solutions with an electric field. The U{1) 
gauge field is simply A = ^dt and we parametrize our metric in four spacetime dimensions 
as: 

= ^^^^dt^ - i id-' + dy') - j^.dz^ (4) 

The radial coordinate is defined for z > 0, where z = is the location of AdS boundary. 
All coordinates are dimensionless, according to (II]). This form of the metric is sufficiently 
general to model any configuration of static and isotropic charged matter. Development of 
a horizon at finite z is signified by the appearance of a zero of the function f{z), fi^n) = 0. 
From now on we will set L = 1. 

We will now proceed to derive the equation of motion for the Dirac field. From ([3]), the 
equation reads: 
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In the metric (jl]) we can always eliminate the spin connection 8| by transforming: 
After decomposing into radial projections defined as: 



*± = -(i±r^)vi>, 



(5) 



(6) 



(7) 
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in a basis where = diag(l, 1, —1, —1), the Dirac equation in matrix form becomes: 

7a, M ^ D M . (8) 



Here the matrix D is the differential operator along the transverse coordinates (x, y) and 
time, which we will specify shortly. 

We will give the solution of the Dirac equation in the cylindrical coordinates, which 
will serve as the input to the calculation of bulk fermion density in WKB approximation. 
Introducing the cylindrical coordinates as {t, x, y, z) h-> (t, p, 0, z) we make the separation 
ansatz: 

V(2,P,0)\ f duj ( F{z)Kiip,(f)) 



(9) 



^M/_(z,p,0)y J 2vr \-G{z)K2{p 
where the F, G are are taken as scalars and the modes Ki^2 are in-plane spinors. The Dirac 
equation then takes the form: 

-d/^M {E{uj,z) + M{z))a,\ ( FK, 

'E{co,z)-M{z))as -d/^/M ] \-GK, 




(10) 



We recognize the matrix at the right hand side as D j \fj . The terms E and M have the 
meaning of local energy and mass terms, respectively: 

ph{z)/2 _ 

E{z) = -^^(c. + e<!>{z)), M{z) = -^==. (11) 
The in-plane operator d acts on each in-plane spinor as: 

d=\ \ (12) 




with d = e*'^(9p + d,j)/p). To maintain the separation of variables in ( JTOj) . we require dKi = 
XiKi, where |Aip corresponds the momentum-squared of the in-plane motion of the particle. 
The physical requirement that this momentum be the same for both radial projections 
translates into the condition IA2I = |Ai|. Consistency of the separation of variables then 
shows us that K2 = crsKi and thus Ai = — A2 = k and the reduced radial equation becomes: 




-k E + M 
M -E k 
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(13) 



with k = hj (let us note that Eq. ( IT3l) is for the pair (F, G), whereas the initial equation 
( ITOj) is written for the bispinor (FK\^ —GK2)). For the WKB calculation of the density, it is 
useful to remind that the wave function in Eq. iQ has two quantum numbers corresponding 
to the motion in the (p, 0) plane: A, / (or equivalently the momenta k^, ky in Cartesian 
coordinates). The radial eigenfunctions in 2;-direction provide a third quantum number n. 



III. EQUATION OF STATE OF THE BULK FERMION MATTER 

In this section we construct the model of the bulk fermions in an improved semiclassical 
approximation. We solve the Dirac equation in the WKB approximation, and the density is 
computed by summing a large number of energy levels. This is in the spirit of Thomas- Fermi 
approximation. However, we sum a finite number of exact quantum-mechanical solutions 
for the wave functions rather than approximating the sum by an integral as implied in the 
semiclassical fluid limit. One of the drawbacks of the Thomas- Fermi fluid limit are sharp 
bounds (i.e., discontinuous first first derivative) of density and pressure profiles along the 
radial direction (see e.g. [lO, 12, 34|). As we have already argued, sharp bounds make 



it hard if not impossible to capture several phenomena. In this respect summing WKB 
wavefunctions goes beyond Thomas-Fermi; it both includes quantum corrections as the 
number of occupied states is countably finite and all collective and individual profiles will 
be continuous without sharp edges. It bridges the gap between the all-classical electron 
star jlO| and single-particle quantum mechanical calculation of Dirac hair jl3 |. 



A. WKB hierarchy and semiclassical calculation of the density 

In the framework of quantum-many-body calculations, the first task is to construct the 
induced charge density n{z). Physically, the origin of the induced charge in our model is 
the pair production in the strong electromagnetic field of the black hole. To remind the 
reader, a (negatively) charged black hole in AdS space is unstable at low temperatures, and 



spontaneously discharges the vacuum [25|]. This means that there will be a non-zero net 
density of electrons n{z). One can calculate n{z) in a Hartree approximation as a density 
of non-interacting electrons, compute the collective effect on other fields by this density 
and iterate. Our novel approach is to use WKB methods to efficiently compute the many 
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wavef unctions enumerated by the quantum numbers (A,/,n). 

The algorithm for the WKB expansion of the wave function for Dirac equation is adopted 
from 32|. Even though every single step is elementary, altogether it seems to be less well 
known than its Schrodinger equivalent. We consider the Dirac equation in the form and 
introduce the usual WKB phase expansion: 



<iJ{z) = e^^o'''y^'>Vn^>^(^z) (14) 

with the spinor part x{^)- The phase y{z) can be expressed as the semiclassical expansion 
in 

y{z) = (y^iiz) + yo{z) + yi{z) + ...)■ (15) 
The equations for the perturbative corrections now follow from fll4til5p : 



Dxo = y-iXo, (16) 
DXi = y-iXi + yoXo + Vfd^Xo, (17) 



DXn = y-lXn + \^dzXn^l + ^1 (18) 



1=0 



Notice in particular that y^i/xo is an eigenvalue/eigenvector of D. In our case the matrix 
D has rank two, so there are two eigenvalues/eigenvectors for y_i/xo- y^i and x^- To find 
the first order correction to the phase of the wave function yo, we multiply f|T7|) from the 
left by the left eigenvalue of the matrix D [D is in general not symmetric, so the right 
and left eigenvalues are different): 



= (19) 

(Xo > Xo 



so we can now construct the usual WKB solution of the form \E'-|- = e'^*/^, where q is the 
WKB momentum and 9± the phase. The term y^ is just the first order correction to 9±. 

Finally, let us recall the applicability criterion of the WKB calculation. It is known that 
WKB approximation fails in the vicinity of turning points. The condition of applicability 



^ From the very beginning we put h = \. However, to elucidate the semiclassical nature of the expansion 
we give it here with explicit h. Dirac equation becomes hy/Jdz^ = D^, where ^ = (5'+, ^-), yielding 
the expansion y{z) = + hyo{z) + h^yi{z) + ...). 
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comes from comparing leading and the next to leading term in the expansion (fT5|) : 

y^^'^ «i. (20) 



In terms of E{z) and M{z) introduced in Eq. ffTTl) it gives at k = 0: 

M{z)d,E{z) - E{z)d,M{z 



E{z){E{z) - M{z)) 

1. WKB wave function 



< 1. (21) 



According to ( !T6|) . the leading effective WKB momentum for the motion in z direction 



9 = \y-i\ is: 



q^{z) = E^{z)-M\z)-~k''{z). (22) 



The wave function in radial direction, ^ = {F, —G), is given by the superposition of two 
linear independent solutions 

^{z) = C+x+{z)e'"^'^ + C^xAz)e-''^^\ (23) 

with the phase determined by 

e{z) = r {q{z') + 59{z'))dz' (24) 



kd,k - qd,q +(e-m) (d,E + d,M) 
59{z) = I ^ — . — Uz. (25) 



26| 



2kq 

The constants C+ and C„ are related by invoking the textbook boundary conditions 
for the behavior of WKB wave function at the boundary of the classically allowed region 
(g^(z) > 0) and the classically forbidden region {q'^{z) < 0). The wavefunction in the 
classically allowed region then reads: 

,,,, , C_ I jEiz) + Miz)sm{9{z) - 5e{z)) 

J E{z) - M{z) sine (z) 



59{z) = ArcSin ^ - ^^^^ (27) 



^Je^{z) - M^{z) 



and C is the only remaining undetermined normalization constant. Integrating the probabil- 
ity density over all coordinates in classically allowed region {zi,Z2) gives the normalization 



condition: 



C j'^ dz^^^^- j pdp J d<l)Cl^^a{z, P, P, 0) = 1. (28) 

The metric factor is gzd{z) = g{z)g^^{z), and a{z) is the conversion factor from ([6]). In 
the left-hand side of the equahty we took into account the normahzation of the continuous 
spectrum in the (p, (p) plane. The integration over (p is trivial. The orthogonality relation 
for Bessel functions (which encapsulates the (p, 0) solution) gives the definition of C|^: 

C2I J{\p)J{\'p)pdp = (29) 

and it allows us to express the normalization constant as: 

~ \ "1/2 



q{z) 

where a factor of 27i comes from the integration over (p and an additional factor of 2 from 
the summation over the full four-component wave function, i.e. bispinor (each spinor gives 
E[z)/q{z) after averaging over the fast oscillating phase 6). This completes the derivation 
of WKB wave function and allows us to compute the density. 

2. WKB density 



As in (32I] we find the total density by summing single-particle wave functions in the 
classically allowed region. The WKB wave function is characterized by the quantum numbers 
(A, /, n) with A being the linear momentum in the x — y plane, / - the orbital momentum in 
the X — y plane and n - the energy level of the central motion in the potential well along 
z direction. The bulk density can be expressed as the sum over the cylindrical shells of 
the bulk Fermi surface. Each shell satisfies the Luttinger theorem in the transverse (x — y) 
direction and so the density carried by each shell nxy{z) can easily be found. We can then 
sum over all shells to arrive at the final answer which reads simply J dznz{z)nxy{z). A 
similar qualitative logic for summing the Luttinger densities in the x — y plane was used 



also in [ij] although the model used in that paper is overall very different (see also [? ]). 

Let us start by noticing that the end points of the classically allowed region determine 
the limits of summation over n and A: q^{(jJni A) > 0. Thus, the density in the WKB region 
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is: 

^(^) = ^E/ E / / dppCl\<if{z,p,(j))\'. (31) 

The limit of the sum over the level number n is determined by the requirement that WKB 
momentum be positive; in other words, we sum over occupied levels inside the potential well 
only. Remember that the bulk fields live at zero temperature, hence there is no Fermi-Dirac 
factor. The sum over the orbital quantum number I extends to infinity as the (x, y) plane is 
homogenous and the orbital number does not couple to the non-trivial dynamics along the 
radial direction. We can now invoke the (local) Bohr-Sommerfeld quantization rule: 

j dz^M^ = NwKBTT (32) 

to estimate the total number Nwkb of radial harmonics in the sum. The expression for 
NwKB in combination with ( l30l) then give: 



1 duj„ 



1/2 



Cn = [^-Q^J ' for q{z) » 6e{z), z ^ 1. (33) 

Now we turn the summation over the quantum number n into the integration over energy 
and obtain for the hulk electron density (here we also performed the integration over p using 
the explicit expression for the wave function ( ITOl) and the normalization condition ( l29l) for 
the Bessel functions): 



2 ^ /•2'^ pyj f{z){&(u,,z)~M'ii^z)) Pfii,, E{UJ,Z) 



a(z)2 f - io Jo Jo 47r2g(u;, X, z) 



1=0 

h2. 



After performing first integral over u and then over A we get : 



^(^) = ^^^3^2 (35) 
with qwKB determined by 

q'wKB = E'{0,z)-M'{z). (36) 



^ The given result for n can be compared to the charge density in the electron star limit given in [3l ■ The 
metric functions used there are related to ours as / fe'^/z^ and g H' 1/ fz^, where our metric functions 
are on the right hand side. Likewise, our definition of qwKB is related to kp of [18[ as qwKB = kpj \ff . 



Now the total bulk charge is expressed in 18| as Q = / dzh^[z) where n^(z)'^ n(z)e''/^. In our conventions 



Q = J dzy/^g^^g^^n = J dzn{z)e^^'^ thus giving the same result as in |18|. 
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Notice that this formula corresponds with common knowledge on the density of electron star 



10|. However, even though the formal expression is the same, the self-consistent solution 



for the metric and gauge field is different because of the quantum correction we introduce 
to pressure. The difference is visualized in Fig. lA where we preview our self-consistent 
backreacted WKB star solutions and compare them to the semi-classical (electron star) 
limit. 



B. Pressure and equation of state in the semiclassical approximation 

Following the logic behind the density calculation, we will now calculate the pressure. It 
will actually prove easier to write the equation of state first and then derive the pressure. 
We can start by computing the energy density of the bulk fermions. By definition, it reads 

/•27r poo rfJ-ioc /-Stt pco rfJ-ioc E( z) 

£{z) = V / d(t) dp duju^\z)^{z) = V / # / dp duu ^ 

Jo Jo Jo \i Jo Jo Jo ^'^ 'i^^l 

(37) 

where E{z) is defined in ( fTTl) and the sum limits are the same as in ( IMIl . Performing the 
integration in a similar fashion as when computing n{z) in (13411351) . we obtain 

£ = le$n + - f^M^ ArcSinhi. (38) 
2 2-' M 

Notice that the first term exactly captures the electrostatic energy while the second is the 

one-loop term that encapsulates the quantum fluctuations. The above result is remarkably 

close to the Hartree vacuum polarization correction as it appears in various model energy 

functionals in literature. 

Now the calculation of pressure needs to be done very carefully in our setup. There are 

two possible approaches: 

1. Fundamentally we can express the pressure from the microscopic fermionic Lagrangian 
([3]). By definition it reads 

n,l,\ n,l 

The explicit calculation is tedious but straightforward. Unlike for density case, the 
final sum is not readily performed to obtain a closed-form expression. Instead, we 
integrate numerically over the energy levels Un to obtain the function p{z). 
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0.0 0.2 0.4 0.6 0.8 1.0 



(C) z (D) 

FIG. 1: WKB bulk density n{z) (Eg. 1351 blue lines) and electron star density (red dashed lines). 
Parameter values (A) {^,e,m) = (1.7,1,0.1), (B) {n,e,m) = (1.7,10,1). The classically allowed 
region lies between the turning points and z**, determined by the condition of vanishing WKB 
momentum {q{z^) = ^(z**) = 0). The parameters for (A) are in the classical (electron star) 
regime, with N^kb ^ 1 when WKB approximation is quite accurate. The plot (B) shows a case of 
small N\\rxB where the WKB approximation becomes inadequate and further quantum corrections 
are likely to be important. (C) Bulk density for a range of values (/i, e, m) = (1.7,1,0.1) (red), 
{fi,e,m) = (1.7,5,0.1) (violet), {fi,e,m) = (1.7,10,1) (green) and {fi,e,m) = (1.7,20,1) (blue). 
For large specific charge of the fermion (and therefore a large number of WKB levels in the bulk) 
the solution is dominated by the classically allowed region and looks similar to the electron star 
limit. For smaller e/m values (and thus fewer WKB levels) the quantum correction in the near- 
boundary region becomes more important and the curves are visibly different from the fluid limit. 
(D) Thermodynamical pressure (Eq. I40p . for the same parameter values as in (C). 

2. In a "near" -classical regime, at large occupation number, thermodynamics ought to 
work, so p = dE/dV. While the volume V is difficult to calculate exactly, we can 
obtain a crude estimate in the following way. At unit AdS radius, the volume equals 
the length i = of the classically allowed interval along z axis, i.e. the in- 
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FIG. 2: Comparison between full quantum pressure (dashed blue lines, Eq. l39p and thermodynamic 
pressure (solid black lines, Eq. I40|) for two sets of parameters: (yu/T, e,m) = (1.7,1,0.1) (A) and 
(/z/T, e, m) = (1.7, 5, 1) (B). For comparison we plot also the fluid pressure^ = en<I>/2 (dashed green 
lines) . Expectedly, for Nyyx b ^ 1 (A) the thermodynamic approximation comes close to the exact 
many-body calculation while for N^kb small the level spacing is large and the thermodynamic 
limit is no longer a good approximation to the sum of the contributions of individual levels. Notice 
that both ways of computing pressure yield similar results for large Nwkb but deviate at smaller 

NwKB- 



terval between the zeros of the WKB momentum qwKB{z) = \ E'^{z) — M'^{z). For 



Nwkb ^ 1, it is natural to suppose (and confirmed by the numerics, see Fig. 1) 
that z** is very close to the horizon, is quite far from the horizon and £ ~ 1 — z*. 
For 2 ~ 2;*, it is natural to assume that the electric potential does not deviate much 
from the linear law: $ ^ yu(l — z), because z^, is not far from the boundary. This 
means that the metric function h{z) can be well approximated by a linear function 
h{z) ~ const. (1 — z). Solving the equation q^^s — -^^(^*) ~ M'^{z^) = 0, we get 



(140|) thermodynamic pressure and denote it by pthd to differentiate from the exact 
quantum expression fl39|) . These expression are also the equations of state of the 
system as they connect the pressure to the density. 

The thermodynamical pressure is more convenient for calculations. In spite of its approxi- 




£ ~ 1 — log 



and (EHD gives: 
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mate nature, (HOj) in particular yields a remarkably accurate result when compared with the 
exact quantum pressure at Nwkb ^1-^ 

Finally, it is illustrative to see how we reproduce the electron star pressure [l^ in the 
limit of large density. For n — j- oo, the first term in S and pthd dominates and we obtain 
from dMD and 

Pes = ^e$n (41) 

as expected for an ideal fluid, which corresponds to the electron star approach. The physical 
interpretation of this result (and of the pressure inside the classically allowed region in 
general) is that of a Fermi gas pressure which, as we know, survives also in the limit of 
classical thermodynamics. The comparison of p, pthd and pes is summarized in Fig. [21 
for high and low number of levels. While the thermodynamic approximation fl40p is good 
when Nwkb ^ 1, for small Nwkb both the fluid limit and the thermodynamic limit (??) 
eventually break down and the contributions of individual levels must be taken into account. 



IV. MAXWELL-DIRAC-EINSTEIN SYSTEM 



We have now arrived at the point where we can solve our model self-consistently with 
the Einstein-Maxwell equations. The solution is numerical, using an iterative procedure to 
converge toward the solution. Only in the IR region it is possible to use a scaling ansatz 
to estimate the scaling behavior of the metric and matter flelds, akin to the procedure used 



m 



2l|. This is the flrst attempt at a self-consistent solution including backreaction on the 



geometry with holographic fermions which goes beyond the fluid picture of 101 ]. 

Our calculation is similar to the one for relativistic ideal fluid (i.e. electron star) approxi- 
mation. Because an ideal fluid is dissipation-less one can construct an action as put forward 



m 



16j and used in 



The Lagrangian of this charged fluid coupled to gravity and 



electromagnetism is 



S 



d X 



(d^^f + p 



(42) 



^ We have now calculated the radial pressure, i.e. the fermionie component of the stress tensor . However, 
due to isotropy of the system, the components and wi ll receive the same contribution from the 



fermions. The same happens in the fluid limit, as shown in {lo| . 
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In other words, the contribution of fermions reduces to the pressure p. While we do not take 
the fluid limit in this paper, one can suspect that in the first approximation the influence 
of the corrections to fluid limit {N^j^^ — )■ oo) is fully encapsulated by the correction to the 
classical (or fluid) pressure we found in f l5Sll5Ul) . The emergent isotropy and its implied ideal 
nature of the fluid at large occupation number should ensure this. 

To construct the backreacted geometry, we therefore "replace" the fermionic terms in the 
exact action ([2]) with our effective ideal fluid model in terms of the density and pressure of 
the bulk fermions. The total effective action is represented as S* = S'^ + Sm + Sf, the sum 
of Einstein, Maxwell and fluid part. The only nonzero component of the gauge field is $ 
and the only non-vanishing derivatives are the radial derivatives dz (the others average out 
to zero for symmetry reasons). The nonzero fermion pressure p is that considered in Sec. 
IIIIBI and there is a nonzero (local) charge density 

f^=qn,/g^ = qn-^. (43) 

The fermion fluid term in the effective action thus becomes 

Sf = -J d'x^g (j> + p) . (44) 

Due to the preserved spherical symmetry we may substitute these simplifications directly in 
the effective action to arrive at: 



(45) 



The only components of the stress tensor the fermion kinetic energy contributes to are the 
diagonal ones; the others vanish due to homogeneity and isotropy in time and in the x — y 
plane. From f H5|) we get the equations for the energy-momentum tensor: 

T: = Tl = -\z'e^ + j> + mn + glp. (47) 

With the metric ansatz (jlj), we can now write down our equations of motion: 

{dze-^"dz^) = -j° (48) 

3/ - zdj - 3 = (49) 

3/ - zdj{z) - Szfdzh - 3 = T:. (50) 
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In this article we shall only be interested in finite temperature solutions. The gravitational 
background is therefore a black hole with an horizon: a single zero in the warp function f{z) 
at a finite value z = zh-^ Physically the inescapability of the black hole horizon immediately 
suggests the following boundary conditions. The black hole horizon should have no hair so 
^{zh) = 0; h{z) which characterizes the ratio of the UV and IR speed of light should be 
finite at the horizon: h{zH) = ho- Note that the effective WKB potential felt by the fermions 
blows up at the horizon and that the fermion wavefunctions therefore manifestly vanish at 
zh- This same phenomenon is noted in the electron star at finite temperature which also 



has an "inner" edge outside the horizon 



12 



3^ 



At AdS infinity the boundary conditions are standard in AdS/CFT: for the gauge field 
lim^_j,o $(2;) = /i fixes the chemical potential at the boundary [zq — 0). We normalize 
liiRz^Q f{z) = 1, \imz^Qh{z) = 0. Again the boundedness of the normalized WKB wave- 
functions uniquely fixes the behavior of the fermions. 

Finally, it remains to define the units used throughout the paper. The natural unit of 
energy and momentum is the chemical potential /i and we will express all quantities in units 
of fi. The two thermodynamic parameters are the chemical potential fi and T. As AdS/CFT 
is built on conformal field theories which have no intrinsic scale, the physics only depends 
on the ratio fi/T. 

Let us conclude with an outline of the numerical algorithm, which is not completely 
trivial. The boundary conditions to be implemented are given at different points: some 
are given at the AdS boundary and some at the horizon. Since the system is nonlinear, it 
is necessary to either linearize the system or to shoot for the correct boundary conditions 
with the full nonlinear system. After experimenting with both, we have decided to iterate 
the full, non-simplified system of equations, integrating from the horizon and shooting for 
the conditions at the boundary. The iterative procedure consists of two steps: we start 
with the non-backreacted AdS-RN geometry and compute the density (semiclassical plus 
the quantum corrections) for the the electron charge equal to e/N (where e is the physical 



^ At zero temperature, when the horizon vanishes due to fermionic backreaction (this ineludes also the case 
of Lifshitz geometry) , the boundary condition for / guarantees also the smoothness of the solution on the 
horizon: dzfizn) = 0. This condition ensures that we pick the correct branch of the solution as there are 
typically two families of functions f{z) that satisfy the equations of motion and the condition f[z) = 0. 
One of them has a vanishing derivative whereas the other has finite derivative as z — > 1. 
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charge and some positive integer), then we solve the system of Einstein-Maxwell equations 
f HHUSOj) . afterwards we increase the fermion charge to 2e/N , calculate the charge density in 
the background (/, /i, $) taken from previous iteration and solve for this density Einstein- 
Maxwell equations (14811501) . We repeat this procedure for charge 3e/N, 4e/N etc. After 
iterations we have arrived at the physical value of the charge e. Then we do more iterations 
with fixed charge e to ensure that the solution has converged, checking that the set of 
functions (/, h, $) does not change from iteration to iteration. In this way we achieve the 
self-consistent numerical solution of the Maxwell-Dirac-Einstein system of equations. The 
integration is always done from the horizon, shooting for the conditions for $ and h at the 
boundary, since it is well known that integrating from the AdS boundary is a risky procedure 
as it is next to impossible to arrive at the correct branch of the solution at the horizon. 



V. PHASES OF HOLOGRAPHIC FERMIONS 



We can now analyze the structure of both the bulk and the field theory side as a function 
of the parameters T/fi, e and m. We first shortly discuss the nature of the bulk solution 
for the geometry and gauge field and notice some qualitative properties. The typical way 
that the solutions to the WKB-Fermi-Einstein system (14811501) look is illustrated in Fig. |3l 



The near-horizon scaling o: 
the light of earlier models 



the metric and gauge field is of Lifshitz type, as expected in 



3 



35] . Notice that we are working at finite temperature and 



thus do not impose the IR boundary conditions for the metric functions which correspond 
to the Lifshitz geometry. Our finding of Lifshitz scaling is purely numerical, with the simple 
boundary conditions discussed above. In the figure, we plot also the electron star solution 
for comparison. One should be careful in comparing the two, however, as the electron 
star corresponds to the limit e — ?■ and thus cannot be compared directly (i.e., for the 
same parameter values) to our WKB star. Our convention is to first define the electron 
star by choosing the total charge density Q and the parameter m = m/en, where k, is the 
gravitational constant whose value is fixed by the normalization of the action (142|) . For the 
WKB star, we impose the same value of Q, while the value of m is found as m = men 
(for WKB star we can control e as and independent parameter). Relative proximity of the 
solutions for large A^ seems to confirm that this is a physically meaningful way of comparing 
the the models. 
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FIG. 3: Profiles of the metric functions f{z) (red) and e^^^^ (violet), the gauge field ^(z) (green), 
density n{z) (blue) and the pressure ^(z) (cyan) at zero temperature, for {fi/T,e,m) = (1.7, 1,0.1) 
(A) and for {p/T, e, m) = (1.7, 10, 0.1) (B). Solid lines are calculated from our model while dashed 
lines are the electron star solution for the same parameter values. For better visibility density 
and pressure are rescaled by a constant factor. Near the boundary we always have h{z) — t- and 
^(z) = fi + 0{z), in accordance with the universal AdS asymptotics of the solution but in the 
interior the solutions start to deviate. Most striking is the absence of sharp classical edges in 
density and pressure. The difference in pressure will turn out to be crucial in moving away from 
the fluid limit. Here we have not shown the solution with N^kb = ^- this case deviates from the 
electron star (N\ykb oo) so strongly that it does not make sense to compare it. Indeed, 4 <C oo! 



A. Thermodynamics 

We can now use these full solutions to determine the macroscopic characteristics of the 
dual strongly coupled fermion system. Let us first derive the free energy of the boundary 
field theory. According to the dictionary, it is equal to the (Euclidean) on-shell action, which 
contains both bulk and boundary components: 

p (^on~shell _|_ (^on— shell (51) 

bulk bixd V ' 

We have already discussed the bulk action in the previous section. We will again approximate 
the fermionic contribution (jH]) by its leading term, the pressure. 

In computing the free energy using AdS / CFT a crucial part is often played by boundary 
terms in the action. It encapsulates the regularizing terms that eliminates z — j- divergences, 
enforces a Dirichlet boundary condition for the gauge field, but it also provides the kernel 
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for the fermionic correlation functions 



0,0 



Sbnd = i r-n^F^^A^' + ) , (52) 

JdAdS \^ / 

with h being the induced metric on the boundary {h = -^{—l/f^z = 0), 1, 1)) and and 
\[^_ are radial projections of the wave function as in Eq. (JTj). By dAdS we have denoted the 
boundary of the AdS space. Let us now briefly show why these boundary fermion terms do 
not contribute to the free energy, but that the leading fermion contribution is the (one-loop) 
effective pressure. Essentially the point is that only normalizable modes of the field are 
occupied and hence they cannot contribute to the boundary action as they die off too fast. 



The Dirac field asymptotics at the boundary are given by 13|: 



J^XT^ 5/2+™ ^ _ _ = 5_z3/2+- + . . . (53) 

+ 2m + 1 ' ^ ^ 

At the same time the electromagnetic boundary term reduces to ^dz^\z=o = —fJ'P, where 
p is the total boundary (not only fermionic) charge density, read off from the subleading 
"response" of the bulk electrostatic potential lim2_j.o ^{z) = fi — pz + . . .. The regularized 
boundary action now reads 



Simd = lim S{zq) + lim / d^x 



(54) 



20^0 ' 20^0 J [2{2m + 1) ' " 2' 

Since m > — 1/2 is the fermionic unitarity bound in AdS/CFT, the first term always vanishes 
in the limit Zq — 0. The total on-shell action, i.e. the free energy is therefore 



dzd^Xy 
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-pp (55) 



B. Constructing the phase diagram: Quantum corrections imply a first order 
thermal phase transition to AdS-RN 

The condensed matter context in which we are discussing AdS / CFT is that of an emergent 
finite density fermionic ground state out of an UV CFT. In the deep UV or at very high 
temperatures T/ p the chemical potential should be negligible and we should recover as the 
preferred groundstate the UV CFT at finite T / p. The gravitational dual of this is the AdS- 
RN black hole. It describes a conformal critical phase with no Fermi surfaces. As we lower 
T//i an instability should set in towards a state with a finite occupation number of fermions. 
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In the probe analysis one indeed finds several normalizable wavefunctions signalling the 
existence of states with distinct occupation numbers. They are the bulk counterpart of the 
existence of non- Fermi-liquid Fermi surfaces js], oi Q, 22 1. A crucial qualitative aspect is 



that due to their fermionic nature the wavefunctions of these normalizable modes can never 
"grow" . From a microscopic point of view it therefore appears that any fermion driven phase 
transition cannot be second order. In the fluid limit, however, the transition was found to 
be third order. There is no conflict because new analytic behavior can emerge in the fluid 
scaling limit where the number of Fermi surfaces is taken to infinity.^ It does mean that 
one has to be quite careful in the fluid limit as for fermions these corrections can change 
macroscopic quantities. For any finite number of Fermi surfaces we should discover a first 
order transition. We did indeed find this earlier in the Dirac hair approximation valid for 

NwB = 



1 13|. With the WKB construction put forward here, we will show that this is 
indeed so for any finite Nwkb- 

Fig. m shows the behavior of the free energy F{T/fi) of the WKB corrected star con- 
struction for different parameters e, m, corresponding to a different number of levels Nwkb 
(which roughly equals the number of Fermi surfaces Nwkb)- In the high temperature phase 
the preferred state with lowest F{T/fi) is that of the pure AdS-RN. Since there are no 
occupied fermionic states it is independent of the fermion charge and mass. In the low 
temperature phase the preferred phase is the WKB star. Where the phase transition occurs, 
one immediately sees the characteristic first order cusp in F{T/fi) whose non-analyticity 
indeed becomes clearer as Nwkb decreases. The panel (B) of the figure makes this clear by 
showing the vicinity of the phase transition. 

The first order nature of the phase transition can in fact be understood analytically with 
this WKB construction. The argument is the same as for the exact fluid limit of the electron 



^ Note that there is a crucial subtlety in the fluid limit in AdS/CFT with a flat Minkowski-space boundary. 
Normally one needs a thermodynamic "fluid" limit to even be able to discuss the notion of a phase 
transition. In global AdS, or conventional Tolman-Oppenheimer-Volkov neutron stars, a bound on the 
number of radial modes, implies a countable number of states. However, this is not so in AdS/CFT with 
a flat Minkowski-space boundary. For each radial mode there is still a formal infinite number of modes 
distinguished by the transverse momentum. The phase transition discussed here is where one considers 
N/Vtransverse oo. It restorcs one's intuition that the emergence of each single Fermi surface dual to 
each single radial mode is associated with a macroscopic phase transition. Wc thank Scan HartnoU for 
emphasizing this. 
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star jl2|. Starting from low temperatures, at the transition point the bulk density n vanishes. 
In the WKB construction that means that the turning points coincide: = z**. Assuming 
that the transition is dominated by the behavior of the fermions and that the contribution on 
the geometry change due to backreaction is small near the critical temperature, the relevant 
part of the free energy of the system is given by Fpermi ~ Jn" P- The first, "fluid limit" term 



in the pressure (140!) is proportional to $?i and it is analyzed in detail in 12|. It yields the 
scaling Fpermi ~ {T — TcY. This indicates a third order transition at the semi-classical level. 
The new, second, quantum term will change this, however. The vanishing of the classically 
allowed region means E ^ M in the whole (narrow) region z^, < z < z^:^. One can thus 
expand E = M + 6z x 6E/6z + . . . and analyze the leading terms in 6z. It is easy to see 



that its expansion starts from a constant: $M^e '^f y + e^E'^ = const. + 0{5z), where 
6z Its integral thus scales as Ff^j^Yiii Sz. Now, for a vanishing bulk charged 

fluid/emerging charged black hole, the principle of detailed balance predicts that the charge 
of the former equals the charge of the latter: n6z = ubh^zh, where the charge densities 
of the bulk and the black hole are n and ubh, respectively, and 6zh is the change in the 
position of the black hole horizon. The crucial insight is that the densities can be assumed 
constant for vanishing 6z and Sz^- We thus find 6z ~ 6zh ^ T — T^. We can now write 



Fpermi = -f'/lrmi + ^Fpermi- Wc kuow that -P/g^^^j ~ (T - T^f [jj], but we have now shown 
that 

AFpermi T - T^. (56) 

At the quantum level the transition is always of first order. The quantum correction is sub- 
leading at general T values, but becomes leading as the phase transition point is approached. 
Finally, we remark that, if one considers the bulk free (or internal) energy J dzS given in 
Eq. ( 1371) using the similar scaling reasoning, one arrives at the same conclusion: F T — Tc- 
This confirms the intuition that the bulk and boundary thermodynamics are equivalent at 
leading order, i.e. the difference Fi,uik ~ F does not contain first-order terms in T — T,, and 
thus does not change the order of the transition. 

The numerics just confirms this analytic prediction of a first order phase transition. The 
field theory interpretation of the discontinuous nature of the transition to a phase with Fermi 
surfaces is simple: fermions do not break any symmetry but the discharge of the black hole 
does signify that the ground state is reconstructed due to the formation of a rigid Fermi 
surface. The only way to reconstruct the ground state without breaking any symmetries 
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FIG. 4: Free energy as a function of temperature F{T). The abrupt change of the derivative 
signifies the first order transition between the finite density phase and the pure black hole (with 
zero bulk fermion density), in line with the analytical prediction of the first order transition from 
the second term in the bulk free energy in Sec. IV B[ We show the calculations for three different 
values {fi,e,m) of the system parameters: (1.7,3,0.1), N\ykb{T = 0) = 40 in blue, (1.7,10,0.1), 
Nwkb{T = 0) = 20 in red and (1.7, 10,0.7), Nwkb{T = 0) = 11 in violet. Notice how the slope 
of F in the low-temperature phase decreases as the number of levels increases: for N\ykb oo 
we reach the electron star limit when the transition becomes continuous. Panel (B) shows the 
vicinity of the critical temperature for three sets of parameter values, to make the cusp in F{T) 
clearly visible. In the high temperature (RN) phase the curves F{T) fall on top of each other as 
one expects for the RN black hole with n = 0. The behavior in the low-temperature phase (with 
non-zero density) is different for the three curves as the value of the charge affects the behavior of 
the bulk fermions. For presentation purposes, the curves have been rescaled to the same transition 
temperature; in general, however, {T/ is not universal and will differ for different corners of the 
parameter space. 



is precisely the first order transition of the density van der Waals liquid-gas type. This is 
the macroscopic counterpart to the probe analysis where the Grassman nature of fermions 
Pauli blocks the growing of mode functions. A van der Waals liquid-gas first order type 



transition is indeed seen in 



13| for the first order transition from N^^b = 1 Dirac hair 



state to AdS-RN. The confusing point was that electron star/AdS-RN transition valid in 



the strict Nv/kb oo fluid limit was found to be third order 
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34j . Here we show that 



this change in the nature of the phase transition is an artifact of this Nwkb oo limit. 
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Instead the expected first order behavior is recovered for any finite value of Nwkb- 



VI. DISCUSSION AND CONCLUSIONS 
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FIG. 5: Applicability of various approximations as a function of the ratio of the fermion charge and 
the total charge of the system, Q/e: Dirac hair, confined Fermi liquid, our present WKB-model, 
electron star. Dirac hair and electron star are the simplest and most flexible approximations but 
limited to the extreme ends of the Q/e axis. 



In this paper we have constructed an improved semiclassical model of holographic 
fermions in AdS4 space, aimed at understanding the phase diagram of strongly coupled 
Fermi and non-Fermi liquids. The model combines a WKB approximation with a Hartree 
summation to approximate a finite N^kb charged fermion state in AdS coupled to both 
gravity and electromagnetism. The dominant effect is a quantum correction to the pressure 
and energy density ("vacuum polarization") of the conventional Nwkb oo classical model 
- the electron star. This finite N-^kb approach has allowed us to address the intermedi- 
ate fermion charges which cannot be modeled satisfyingly with any of the previously used 
models. 

By studying the free energy of the system we can now construct the full phase diagram of 
the system. Most importantly, we find a universal first order phase transition from a finite 
density to a zero density (Reissner-Nordstrom, quantum critical) phase. The discontinuity of 
the density comes from the quantum term in the internal energy. This term is always present 
but its relative contribution to the free energy decreases with the inverse of the number of 



radial modes Nwkb- The extreme limit N\ 
continuous phase transition found in 12 



w 



KB ^ oo reproduces the unexpected third order 



34j . Nevertheless, in any real system with finite 



fermion charge the discontinuity will be present, which fits into the general expectation 
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that the thermal phase transition of a fermionic system should be of the van der Waals 
(liquid-gas, Ising) type. 



How to relate the electron star 



10[ |. Dirac hair [13| and the (confined) Hartree Fermi 



liquid [l^ to our new phase diagram? All models use the same microscopic action for a 
Dirac fermion with charge e and mass m, but the system is approximated in different ways. 
The electron star is the fluid limit of the equations of motion, yielding the Openheimer- 



Volkov equations in the bulk. As explained in 



19|, this approximation is valid in the limit 



of infinite occupation number Nwkb — ?■ oo, e — > with the total charge density fixed 



Q = NivKB^- In addition, the mass m — )■ while m = m/ y/WwKBe is fixed. The Dirac hair 
departs from the opposite limit, treating the bulk fermion as a single collective excitation 
with Nwkb = 1. To obtain a macroscopic charge density one essentiall y h as to take e ^ 1. 



Finally, the confined Fermi liquid of 



14| and its deconfined version 17| improve on the 



Dirac Hair by using a standard Hartree summation of the non-interacting bulk Fermi gas. 
It works for all Nwkb ~ 0(1) and this significantly increases the region of applicability but 
at the cost of substantial practical complications, in particular if one wishes to take into 
account the backreaction on the metric Our model takes a similar Hartree approach 

but simplifies the wavefunction calculation by using the WKB approximation. This inherently 
assumes semiclassical dynamics and large number of energy levels Nwkb ^ 1 in the bulk. 
Nevertheless, we do not make the assumption of zero energy spacing Nwkb — ^ oo necessary 
for the fluid approximation: our model thus works well in the intermediate regime where 
Nwkb is finite but large compared to unity. This message is illustrated in Fig. [5l emphasizing 
the singular nature of both the electron star and the Dirac hair. 

The natural next step departing from this WKB treatment is to employ a fully quantum- 
mechanical density functional method. It is, in fact, not a significant complication compared 
to the approach of this paper: the recipe for computing the density n will be replaced by 
a somewhat more complicated functional of the gauge field and the metric, which needs 
to be computed iteratively. We anticipate that this will not alter the qualitative picture, 
although the number might change significantly. The main conclusion of our paper is that the 
singular fluid limit of bulk fermions when coupled to AdS gravity can lead to macroscopically 
anomalous results. Finite Nwkb corrections are crucial to get the correct answer. 
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